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Abstract. We study Lorentzian wormholes in the ghost-free bigravity theory de¬ 
scribed by two metrics, g and f. Wormholes can exist if only the null energy condition 
is violated, which happens naturally in the bigravity theory since the graviton energy- 
momentum tensors do not apriori fulhll any energy conditions. As a result, the field 
equations admit solutions describing wormholes whose throat size is typically of the 
order of the inverse graviton mass. Hence, they are as large as the universe, so that in 
principle we might all live in a giant wormhole. The wormholes can be of two different 
types that we call W1 and W2. The W1 wormholes interpolate between the AdS spaces 
and have Killing horizons shielding the throat. The Fierz-Pauli graviton mass for these 
solutions becomes imaginary in the AdS zone, hence the gravitons behave as tachyons, 
but since the Breitenlohner-Freedman bound is fulhlled, there should be no tachyon 
instability. For the W2 wormholes the g-geometry is globally regular and in the far 
held zone it becomes the AdS up to subleading terms, its throat can be traversed by 
timelike geodesics, while the f-geometry has a completely different structure and is not 
geodesically complete. There is no evidence of tachyons for these solutions, although a 
detailed stability analysis remains an open issue. It is possible that the solutions may 
admit a holographic interpretation. 
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1 Introduction 


Lorentzian wormholes are hypothetical held-theory objects describing bridges connect¬ 
ing different universes or different parts of the same universe. They could supposedly 
be used for momentary displacements over large distances in space. In the simplest 
case, a wormhole can be described by a static, spherically symmetric line element 

ds^ = —Q‘^{r)dt‘^ + dr'^ + R‘^{r){d'd‘^ + sin^ 'ddif'^), (1.1) 


where Q{r) = Q{—r) and R{r) = R{—r), both Q and R are positive, and R attains a 
non-zero global minimum at r = 0. If both Q and R/r approach unity at inhnity, then 
the metric describes two asymptotically flat regions connected by a throat of radius 
1?(0). Using the Einstein equations Gj) = 87rGT^, one hnds that the energy density 
p = —Tq and the radial pressure p = in the throat at r = 0 satisfy 


R" 


< 0 , 


P = 


< 0 . 


P + P 


AtiGR 


SnGR^ 


( 1 . 2 ) 






It follows that for the wormhole to be a solution of the Einstein equations, the matter 
should violate the null energy condition {T^yV^v'^ > 0 for any null v^). This shows that 
wormholes cannot exist in ordinary physical situations where the energy conditions are 
fuiailed. 

However, it was emphasized [1, 2] that wormholes could in principle be created 
by the vacuum polarization effects, since the vacuum energy can be negative. This ob¬ 
servation triggered a raise of activity (see [3] for a review), even though the wormholes 
supported by the vacuum effects are typically very small [4, 5] and cannot be used 
for space travels. Another possibility to get wormholes is to consider exotic matter 
types, as for example phantom fields with a negative kinetic energy [6-8]. Otherwise, 
one can search for wormholes in the alternative theories of gravity, as for example in 
the theories with higher derivatives [9, 10], in the Gauss-Bonnet theory [11-14], in 
the brainworld models [15], or in the Horndeski-type theories [16] with non-minimally 
coupled helds [17-20]. 

A particular case of the Horndeski theory is the Galileon model [21], which can be 
viewed as a special limit of the ghost-free massive gravity theory [22] . This latter theory 
has recently attracted a lot of attention (see [23, 24] for a review), because it avoids 
the long standing problem of the ghost [25] and could in principle be used to describe 
the cosmology. In particular, it admits self-accelerating cosmological solutions and 
also black holes (see [26, 27] for a review). At the same time, wormholes with massive 
gravitons have never been considered. To hll this gap, we shall study below wormholes 
within the ghost-free bigravity theory. 

The ghost-free bigravity [28] is the extension of the massive gravity theory con¬ 
taining two dynamical metrics, g^y and f^y. They describe two gravitons, one massive 
and one massless, and satisfy two coupled sets of Einstein’s equations, 

G^y{g) = K, mg, /), G^yif) = Uig, /), (1.3) 

where and are the graviton energy-momentum tensors. What is interesting, 
these tensors do not apriori fulhl the null energy condition [29], which suggests looking 
for wormhole solutions^. 

There are two possible ways to interpret the two metrics in the theory. One 
possibility is to view them as describing two geometries on the spacetime manifold. 
Each geometry has its own geodesic structure, and in principle one could introduce two 
different matter types - a g-matter that follows the g-geodesics and does not directly 
see the f-metric, and an f-matter moving along the f-geodesics. However, it is not 
always possible to put two different geometries on the same manifold, and in fact we 
shall present below solutions for which the spacetime manifold is geodesically complete 
in one geometry but is incomplete in the other. We shall therefore adopt the viewpoint 
according to which only the gfj^y describes the spacetime geometry, while the f^y is a 
spin-2 tensor held whose geometric interpretation is possible but not necessary. 

^Even though the energy conditions are not fnliilled, this does not necessarily mean that the energy 
is negative, and in fact the analysis in the massive gravity limit indicates that the energy is positive 
in the asymptotically flat case [30, 31]. 
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In what follows we shall study the wormhole solutions in the system (1.3). It turns 
out that such wormholes exist and are gigantic, with the throat size of the order of 
the inverse graviton mass, which is as large as the universe. Therefore, if the bigravity 
theory indeed describes the nature, we might in principle all live in a giant wormhole. 
We hnd wormholes of two different types that we call W1 and W2. For the W1 solutions 
both metrics interpolate between the AdS spaces, and the g-geometry is either globally 
regular (the Wla subcase) or it exhibits Killing horizons shielding the throat (the Wlb 
subcase). The Fierz-Pauli graviton mass computed in the AdS far held zone turns out 
to be imaginary, so that the gravitons behave as tachyons. For the Wla solutions the 
graviton mass violates the Breitenlohner-Freedman (BF) bound, hence these solutions 
must be unstable, but for the Wlb solutions the bound is fulhlled, which suggests 
that the tachyon instability is absent. For the W2 solutions the g-metric is globally 
regular and in the far held zone becomes the AdS up to subleading terms, while the f- 
geonietry has a completely diherent structure and is not geodesically complete. There 
is no evidence of tachyons for these solutions, although a detailed stability analysis 
remains an open issue in all cases. It is possible that the solutions may admit a 
holographic interpretation. 

The rest of the paper is organized as follows. In section 2 we introduce the ghost- 
free bigravity theory, whose reduction to the spherically symmetric sector is given in 
section 3. The master held equations and their simplest solutions are presented in 
sections 4 and 5. The local solutions in the wormhole throat are obtained in section 
6, while section 7 presents the global solutions. The geometry of the solutions, their 
global structure, geodesics, etc, are considered in section 8. Other properties of the 
solutions, in particular their stability, are briehy discussed in the hnal section 9. 

2 The ghost-free bigravity 

The theory is dehned on a four-dimensional spacetime manifold endowed with two 
Lorentzian metrics gn„ and with the signature (—, -f, -f, -f-). The action is [28] 

^ = ('-I* 

where m is the graviton mass and the two gravitational couplings fulhll ki + K 2 = 
The metrics and all coordinates are assumed to be dimensionless, with the length 
scale being the inverse graviton mass 1/m. The interaction between the two metrics is 
expressed by a scalar function of the tensor 

7 ^ = ( 2 . 2 ) 

where is the inverse of gf^i, and the square root is understood in the matrix sense, 

i.e. 

= g^'^fau- (2.3) 

If Aa (A = 0,1, 2, 3) are the eigenvalues of then the interaction potential is 

4 

U = (2.4) 

n=0 
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(2.5) 


where bh are parameters, while ttkij) are dehned by the relations 
Wo(7) = 1, 

^1(7) = = [7], 

A 

1 / 2 ( 7 ) = E - 71), 

A<B 

^ 3 ( 7 ) = Y1 = ^([ 7 ]^ - 3 [ 7 ][ 7 l + 2 [ 7 ^]), 

A<B<C 

1 / 4 ( 7 ) = A„A,A2A3 = i(|7|' - 6|7|"|7"I + 8|7||7’] + 3(7? - 6 I 7 I). 

Here, using the hat to denote matrices, one has dehned [7] = tr(7) = 7^^, [7^] = 
tr(7^) = (7^)'^^. The two metrics actually enter the action in a completely symmetric 
way, since the action is invariant under 

g^^l' ^ fiMu, ^1 7^ ^ 2 , bk -H- b4_k ■ (2.6) 

The parameters bk can apriori be arbitrary, but if one requires the hat space to be a 
solution of the theory and m to be the Fierz-Pauli mass of the gravitons in hat space, 
then the hve bk are expressed in terms of two arbitrary parameters, sometimes called 
C3 and C4, as 

^o = 4 c 3 + C4 —6, 6i = 3 —803 —C4, 62 = 203 + 04 — !, 63 = ~(c3 + C4), 64 = 04. (2.7) 

Varying the action (2.1) with respect to the two metrics gives two sets of the Einstein 
equations (1.3), where and T)) are obtained by varying the interaction, 

T't= rC = -^T", (2.8) 

with 

r't = {61 Wo + 62 Wi + 63 W2 + 64 W 3 } 7 't 

— {62 Wo + 63W1 + 64 W 2 }( 7 ^)^j, 

+ {63 Wo + 64 Wi}(7^)'^j, 

-&4Wo(7T.. (2.9) 

The Bianchi identities for the Einstein equations imply that 

( 9 ) „ if) 

Vprf = o, Vpr/ = o, (2.10) 

( 9 ) if) 

where Vp and Vp are the covariant derivatives with respect to g^y and f^y. 

We consider the vacuum theory, but one could also add a matter. The equivalence 
principle and the absence of the ghost require the matter to be coupled to one of the 
two metrics but not to both of them at the same time. One could also introduce a 
g-niatter for the g-metric and an f-matter for the f-metric. This is important in what 
follows: test g-particles will follow geodesics of the g-metric and will not directly feel 
the f-metric. 
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3 Spherical symmetry 


Let us choose both metrics to be spherically symmetric, 

dsl = —Q‘^dt^ + + R^dVP' , 

y /\ z 

dsj = -q^df + ^ + U'^dVL^, (3.1) 

where Q, A, i?, g, PL, U depend on the radial coordinate r and dVl? = d'd'^ + sin^ 'ddip^ . 
In this case the tensor 7^^, in (2.2) becomes 


7 


V 


diag 


^ ^ U U 

Q'W'R'R 


The formulas (2.8) then give 


T^, 


7-M 

' U 


diag[T»,ri',T|,r|], 
diag [7.7,7.7?], 


where 


rpO 

-‘-0 

T/ 

rTi2 

-^2 

u7f 


Here u = U/R and 


-Vo-Vi—, 

-Vo-VA, 

-n - -Pi 

Q 

-V2 - Vi 

Q 


-Vo 


qA 

'W 


-Vo — Vo 


q A 


Vi 


QW 

qA 


Rm — bm T 2,bjjij^iU T byyi^2U , 

Rm T ^m+1^ (^ 0, 1, 2). 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


As one can see, the energy-momentum tensors do not apriori fulhll any positivity 
conditions. The independent field equations are 

G^oig) = ^iT^, 

G\ig) = K^Tl, 

G°(/) = 72'C, 

Gl{f) = K2T,\ (3.6) 
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plus the conservation condition V/i = 0 , which has only one non-trivial component, 
V„n‘=(r;)'+ +2 =0. (3J) 

The conservation condition for the second energy-momentum tensor also has only one 
non-trivial component, 

V,. t;" = + ^ (- ro") + 2 ^ y - 7? j = o, (s^s) 

but this condition is not independent and actually follows from (3.7). As a result, 
there are 5 independent equations (3.6),(3.7), which is enough to determine the 6 field 
amplitudes Q, A, R, q, W, U, because the freedom of reparametrizations of the radial 
coordinate r —)■ r(r) allows one to fix one of the amplitudes. 


4 Field equations 

Let us introduce new functions 


N = AR', Y = WU ', 
in terms of which the two metrics read 


(4.1) 


dsl = -Q^dR + 


dR^ 

'W 


+ R^dn^ , 


dsj = —q^dR + -I- U'^dVR. 


(4.2) 


The advantage of this parametrization is that the second derivatives disappear from 
the Einstein tensor, and the four Einstein equations (3.6) become 


r' = {R'YVi + U'NV2) 


2 UNY 

Q' = ~ (ki{QVo + qRi) + 


2UY 
Q{N^ -1)\ RR' 


q' = - {^2 {QVi + qV2) + ^2 


i?2 J 2A2 ’ 
g(y 2 -l)^^ R^U' 
2Y^U ■ 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


The conservation condition (3.7) reads 


V,TR{g) = ^ 


1 - 


N 

Y 


q 


dVo + — dVi ) -|- 


q' NQ'U' 


Q YQR' 


- n = 0, (4.7) 
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and using Eqs.(4.5),(4.6), this reduces to 


where 


R^QV„Tng) = yC^o, 


(4.8) 


c 


/ B'Pf 
V ^ 2UY 


- 1^1 


R^VoVi 

2N 


f R^ViV 2 R^Vf 

I ^2-— ^1-H“ 

V 2UY 2N 


{N^ - 1) RVi 
2N 

(F 2 _ 1 ) 

2UY 


+ {N- Y)RdVoj Q 
+ {N-Y)RdVi] q, 


(4.9) 


with 

dPm = 2 {bm+i + bm+2u) (m = 0,l). (4.10) 

The conservation condition (3.8) becomes 

(f) 

-U\v,Tnf) = j;^C = 0. (4.11) 

The two conditions (4.8) and (4.11) together require that either U' = R' = 0, in which 
case both metrics are degenerate, or that 


C = 0. (4.12) 

As a result, we obtain the four differential equations (4.3)-(4.6) plus the algebraic 
constraint (4.12). The same equations can be obtained by inserting the metrics (4.2) 
directly to the action (2.1), which gives 

/\/f2 r 

S = An —/ Ldtdr, (4-13) 

J 

where, dropping the total derivative, 

1 f{l-N^)R' 

V N 

QR^R' fQR" 


2RN']q + 

J K2 

U' qR^R'\^ 

+ \Vi 


1 [{l-Y^)U' 


N 


Y 

qR^U' 

Y 


-2UY' 


Vo 


(4.14) 


Varying L with respect to iV, V, Q, q gives Eqs.(4.3)-(4.6), while varying it with respect 
to R,U reproduces conditions (4.8) and (4.11). The equations and the Lagrangian L 
are invariant under the interchange symmetry (2.6), which now reads 


K2i Q q, N V, R f/, byn -f-)- ^4—m • (4-15) 

Equation (4.3)-(4.6) contain U\ but so far the expression for U' is missing. To obtain 
it, the only way is to differentiate the constraint, which gives 


dC ,,, dCdC dC , dC dC 

5)7^ - 0 . 


(4.16) 


7 















Since the derivatives N', Y', Q', q' expressed by Eqs.(4.3)-(4.6) are linear fnnctions of 
U', this gives a linear in U' relation, which can be resolved to yield 

U' = VuiN,Y,Q,q,R,R\U). (4.17) 

This equation and Eqs.(4.3)-(4.6) comprise together a closed system of hve differential 
equations for hve variables N, Y, Q, q, U. The i?-amphtude is determined by hxing the 
gauge, for example R = r or R' = N. One can integrate the hve diherential equations 
by imposing the constraint C = 0 only on the initial values, and then it will be fulhlled 
everywhere. 

Alternatively, one can integrate only the four equations (4.3)-(4.6) assuming that 
U' in their right hand side is given by (4.17), while U is obtained by resolving the 
constraint. 

Yet one more possibility is to use the fact that the constraint is linear in Q,q. 
Therefore, it can be resolved with respect to q, 

q = EiN,Y,R,U)Q. (4.18) 

Injecting this to Eqs.(4.3),(4.4),(4.17) gives a closed system of three diherential equa¬ 
tions 


N' = Vn{N,Y,U,R,R'), 

Y' = Vy{N,Y,U,R,R!), 

u' = Vu{N, y, U, R, R'), (4.19) 

and when their solution is known, the amplitude Q is obtained from equation (4.5) 
which assumes the form 

Q' = FQ (4.20) 

with 

F = -(..(P. + EP,) + Y^)^. (4.21) 


5 Simplest solutions 

Some simple solutions of the held equations can be obtained analytically [32],[33]. 
They can be of two diherent types described below in this section. They are not of 
the wormhole type, but the wormholes constructed in the next sections approach these 
solutions in the far held zone. 

5.1 Proportional backgrounds 

Let us choose the two metrics to be conformally related [32],[33], 

dsj = , (5.1) 

with a constant A. This implies that 


g = AQ, U = XR, 


Y = N. 


( 6 , 2 ) 




This also implies that Vm = 'Pm{^) are constant. Imposing the Schwarzschild gange, 
R' = 1, the held eqnations (4.3)-(4.6) and the constraint (4.12) reduce to 

[RN^y = l-K^{Vo + \Vi)R\ =0, (5.3) 

and to the condition for A, 

«i(Po + APi) = y (Pi + XV 2 ) = A(A). (5.4) 

This is an algebraic equation which can have up to four real roots. Choosing a root A, 
the solution of (5.3) is 

N'^ = 1- - -^ , Q = const. X N, (5.5) 

R 3 

where M is an integration constant. Depending on value of A(A), this corresponds 
either to the Schwarzschild or to Schwarzschild-(anti)-de Sitter geometry. If the pa¬ 
rameters fefc are chosen according to (2.7), then the equation (5.4) always has a root 
A = 1 , in which case A(A) = 0. 


5.2 Deformed AdS 


Let us set in the equations U' = q' = [32]. This solves Eqs.(4.6),(4.8), while Eqs.(4.3)- 

(4.5) reduce (with R' = 1) to 


{RN‘^y 



Y' 


1 — kiR Vo , 
RVi 


2UN 


Pi, 


(5.6) 


whose solution is 


Q 

Y 


2 M 

- —- KibiUR 

it 

^dR + AN, 



(5.7) 


where M, A, Yo are integration constants. Interestingly, these expressions can describe 
a wormhole geometry, because if 60 < 0 then —)■ -|-cx) for R —)■ cx), while the constant 
M can be chosen such that vanishes at R = h and A^ > 0 for R> h. Introducing 
the radial coordinate 


dR 

VCR) 


r 


■R 


(5.8) 











and setting in (5.7) A = 0, the g-nietric becomes 



geometry. Unfortunately, Eq.(5.7) does not describe an exact solution, because the 
constraint C = 0 is not fulhlled and so the conservation condition (4.11) for the f- 
metric is not satished. However, the leading terms in Eq.(5.7) describe the asymptotic 
form of a more general solution whose amplitudes U, q are not identically constant but 
approach constant values at large R. Specihcally, expanding the held equations at 
large i?, one hnds the following asymptotic solution. 


jV" = -K,h - Ki6iC7^fl + 0(1) = Nl + 0(1). 



Q = ^fl + 0(1)=Q„ + 0(1), 



(5.9) 


with constant f/oo, qoo- Comparing with ~ ^ ~ AR^/S where A = —nibo < 0, 
one can see that the g-metric is the AdS in the leading order, but the subleading terms 
do not have the AdS structure. 

We shall see below that the wormholes approach for R ^ oo either the propor¬ 
tional AdS solutions (5.2),(5.4),(5.5) or the deformed AdS solutions (5.9). We shall 
call these wormholes, respectively, type W1 and type W2. 

6 Wormholes ~ local behavior 

Since we are unable to obtain the wormhole solutions analytically, we resort to the 
numerical analysis. As a hrst step, we impose the reflection symmetry. Let us return 
for a moment to the parametrization (3.1) and require the two metrics to be symmetric 
under r —r. 


Q{r) = Q{-r), A(r) = A(-r), R{r) = R{-r), 
q{r) = g(-r), W{r) = W(-r), U{r) = U{-r). 


( 6 , 1 ) 


Passing then to the parametrization (4.1), it follows that the functions N,Y dehned 
by (4.1) should be antisymmetric. 


N{r) = -N{-r), Y{r) = -Y{-r). 


( 6 . 2 ) 
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This suggests a local power-series solution around r = 0, 


N = Nir + + ... Q = Qo + Q2r^ + ■■■ R = h + R2r^ + ... 

Y = Yir + Y^r^ + ... q = qq + q^r'^ + ... U = ah + U 2 r^ + ... . (6.3) 

Here h = R{0) is the radius of the wormhole throat measured by the first metric, and 
a = U{0)/R{0) is the ratio of the throat radii measured by the two metrics. 

From now on we shall adopt the gauge condition 

N = R', (6.4) 


which implies that 

ds] = + dr^ + R^dn ^, (6.5) 

so that r is the proper distance measured by the g-metric. The next step is to impose 
the held equations to determine the coefficients in (6.3). To begin with, one notices 
that when inserting (6.3) to Eqs.(4.3)-(4.6), (4.12) and expanding the result over r, 
the leading terms are given by Eqs.(4.5),(4.6) which contain a pole in the right hand 
side due to the terms R'/N"^ ~ 1/r and U'/Y'^ ~ 1/r. For the equations to be fulhlled 
in the leading order, the coefficient in front of the pole should vanish, which imposes 
the conditions 

0 — Qo + KiP 1 go = 0 , 

^^2P2 ~ % + ^2P1 Qo = 0) (6.6) 


with Pm = Vm{cr) = bm + ‘2bm+iO' -|- bm +2 ■ Thesc two linear equations will have a 
non-trivial solution if only their determinant vanishes. Therefore, one requires that 

(ki/i^Po — l) [K2h‘^P2 ~ l) ~ KiK2h'^Fl = 0. (6.7) 

If this condition is fulhlled then the solution of ( 6 . 6 ) is 

go = O'Qo 

with 

1 — Kih‘^Po K2h‘^Pl 
Kih'^Pl 1 — K2h‘^^2 
The value of Qo is irrelevant, as it can be changed by rescaling the time, hence we 
set Qo = 1. Eq.(6.7) plays the key role in our analysis and provides the necessary 
condition for the wormholes to exist. For a given h, this is a fourth order algebraic 
equation for a. If we assume for a moment that the parameters bk are given by (2.7) 
with C 3 = C 4 = 0 , then the equation becomes quadratic and gives 

3Kih{K2h‘^ - 1) ± \/i^i{.i^2h‘^ + l)(3h2 - 1) 

- 1 ) -^ 


( 6 . 8 ) 

(6.9) 


- 11 - 






This expression, assuming both ki and K 2 to be positive, will be real-valued if the 
square root is real, which requires that h > l/\/3. Since h is measured in unites of 
1/m, which is of the order of the Hubble radius, it follows that the wormholes are 
gigantic, as large as the universe. 

Let us return to the expansion of the equations over r. Having removed the r~^ 
terms in Eqs.(4.5),(4.6), the next-to-leading order terms are provided by Eqs.(4.3),(4.4), 
which reduce in the order to 


Ni 




-y/^(Po + 


2 C 2 \ 1 


^2 h , ^ 


2 U 2 

Yi 


+ 


Uo 


ahYi 


( 6 . 11 ) 


These relations can be used to express R 2 = Ni/2 and U 2 in terms of hi, while 
Eqs.(4.5),(4.6) considered in the Y order provide similar expressions for Q 2 and 52 - 
Altogether this gives 


Ni 

U 2 

Q2 

(I 2 


2 R, = -^Y„ 


? I 1 + 


«2 

2 aYi 


K2hPi 
Kl (2U2 


Yi, 




Hi 

Y, 


2 Q 2 + 


2R2 - El 

hPi 


(dPo -l- adPi)Qo ) , 


( 6 . 12 ) 


whereas the value of Yi is hxed by the constraint (4.12), 


K 2 h‘^Pi{K 2 + Kia^){dPo -|- 2adPi + a‘^dP 2 ) — 2aPi{K2 -l- kio^) 
2 ah [(^2 + 2 Kiaa)dPo + 2 Kiaa‘^dPi — K 2 a‘^dP 2 ] — 2 ahPi{K 2 + kiCt^) 


(6.13) 


Continuing this process would allow one to recurrently determine all higher order co¬ 
efficients in the expansions (6.3). For example, in the next two orders Eqs.(4.3)-(4.6) 
determine R 4 , Q 4 , Q 4 , in terms of Y^, while the latter is determined by the con¬ 
straint. For given values of the couplings ki, K 2 and bk, the only free parameter in the 
expansions is the wormhole radius h, all other coefficients being hxed by the equations. 
Therefore, the wormholes are characterized by only one continuous parameter, their 
size h. However, since the algebraic equation (6.7) can have several roots, there could 
be several different wormholes with the same radius h. 

It is worth noting that the expressions in (6.12),(6.13) still exhibit the interchange 
symmetry (4.15), even though this is not completely obvious now, when the gauge is 
hxed. Indeed, our gauge choice is grr = 1, whereas directly interchanging the metrics 
would give the solutions in the diherent gauge, where frr = 1- Let us introduce the 
radial coordinate 

r U' 

z = j ydr, (6.14) 
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so that the derivative of a function / with respect to ^ at r = ^ = 0 is 


f = —f ■ 

Jz 2Uo-’^ 


If one interchanges the two metrics, one will have ■H- /^. Since N Y, R ■H- U and 

Q ■H- q, it follows that the coefficients in (6.12),(6.13) should fulfill the relations 


u 

2U2 

Xt 

AUi 


Ni 

R 2 


V, 

U2, 


V 

2U2 

xx 

4f/| 


u = 


U 2 


iVi, 

R 2 , 


(6.15) 


and similarly for Q 2 and q 2 . Here the underlined expressions should be evaluated 
for the interchanged parameter values: ki ■h- K 2 and bk 0 b^-k- A straightforward 
verification shows that the expressions (6.12),(6.13) indeed fulfill the relations (6.15). 


7 Wormholes global solutions 

Skipping the important issue of convergence of the power series in (6.3), the above 
results indicate that the wormhole solutions exist at least locally, in the throat. The 
next step is to construct them globally. To this end, we extend the local solution 
(6.3) towards large values of r numerically, using the standard integration procedure 
described in [34]. Our results are as follows. 




Figure 1. Profiles of the type Wla wormhole solution obtained for the parameter values 
(7.1). The amplitude Q is everywhere positive but q changes sign. For large r the solution 
approaches the proportional AdS background. 

Choosing some values for the parameters ki, K 2 and bk, it turns out that the local 
solution (6.3) extends to the whole interval r G [0, cxo) only for narrow sets of values 
of h. These latter are selected by the condition that a determined by (6.7) is real. In 
addition, the second derivative i?"(0) = 2 R 2 should be positive, since otherwise R{r) 
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vanishes at a finite valne of r. Finally, even if these two conditions are fnlfilled, the 
solntion may exhibit a singnlarity at a finite proper distance away from the throat at 
a point where the derivatives U' or Y' diverge. However, all these problems can be 
avoided (for some parameter valnes) by adjnsting the throat radins h. 

For properly chosen valnes of h the solntion extends np to large r and approaches 
for r —)■ cxo either the proportional AdS backgronnd or the deformed AdS backgronnd 
described in Section 5. According to their asymptotic behavior, we shall call these 
wormholes, respectively, either type W1 or type W2. 

7.1 Type W1 wormholes 

For these solntions the two metrics in the far field zone become proportional to each 
other and approach the proportional AdS backgronnd described in Section 5. However, 
before reaching this asymptotic, either Q or g or both change sign. If only one of these 
amplitndes vanishes, then we say that the solntion is of type Wla. If both Q and q 
flip sign, then the solntion is called type Wlb. 




Figure 2. Profiles of the Wlb wormhole solution obtained for the parameter values (7.4). 
For large r the solution approaches the proportional AdS background. Both Q and q change 
sign. 


In Fig.l we present an example of the Wla solntion obtained by choosing 

Ki = 0.688, K2 = 0.312, C 3 = 3, 04 = — 6 , h = 2.20, a = 0.444, (7.1) 

and with = 6 ^( 03 , 04 ) given by (2.7). The g-metric shows the wormhole throat and 
is globally regnlar, for large r the whole solntion approaching the proportional AdS 
backgronnd with = X^g^y. To see this, we notice that assnming the valnes (7.1), 
Eq.(5.4) gives three possible options for the proportionality parameter A, 

Ai = 1 , A 2 = -1.264, A 3 = 0.358, (7.2) 

with the corresponding valnes of the cosmological constant 

A(Ai) = 0, A(A2) = -7.474, A(A3) = -0.170. (7.3) 
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The numerical solution chooses the last of these three options as the asymptotic. In¬ 
deed, the ratio U/R shown in Fig.l approaches at large r precisely the value A 3 , as 
does the ratio q/Q, while the ratio Y/N (with N = R') approaches the unit value. All 
this agrees with Eq.(5.2). Next, according to (5.5), the amplitude should approach 
the AdS value Nq = 1 — R{\^)R?/?>, and indeed the ratio N/Nq approaches unity, as 
seen in Fig.l. Finally, the ratios Q/R and q/R should approach constant values, which 
is indeed the case. 

The amplitude Q is everywhere positive, but q changes sign at some point, so 
that the metric coefficient /oo = —q‘^ develops a double zero. This corresponds to a 
Killing horizon of the f-geometry and, as we shall see below, the curvature diverges at 
the horizon. At the same time, nothing special happens to the g-metric at the point 
where q vanishes. The g-geometry is everywhere regular and interpolates between two 
AdS asymptotics as r varies from —00 to -|-oo, passing at r = 0 through the wormhole 
throat of size h. The test particles of a g-matter coupled to the g-metric will therefore 
see a regular wormhole. 

It is possible that for some special parameter values there could be solutions for 
which both Q and q are sign dehnite, but we could not hnd them. On the contrary, 
we hnd solutions for which both Q and q change sign, so that both g and f geometries 
exhibit Killing horizons. The g-horizons and the f-horizons are generically located at 
different points, because they are singular, for if they were regular they would coincide 
to each other [35]. An example of the Wlb solution obtained for the parameter values 

Ki = 4.446, K2 = -3.446, C 3 = 1, C 4 = 0, h = 0.426, a = 1.6 (7.4) 

is shown in Fig.2. Notice that ^2 < 0 in this case. Since their both metrics show singu¬ 
lar horizons, one could think that the Wlb solutions are less interesting as compared 
to the Wla ones. However, we shall see below that the Wla solutions are prone to the 
tachyon instability, whereas the Wlb solutions seem to be free of this problem. 

7.2 Type W2 wormholes 

For these solutions both metrics are globally regular and the g-geometry describes a 
wormhole, but the f-geometry is completely different. An example of such a solution 
is shown in Fig.3 for the parameter values 

Ki = 0.574, = 0.425, C 3 = 0.1, C 4 = 0.3, h = 3.731, a = 0.55. (7.5) 

For these solutions the amplitudes U and q approach asymptotically constant values 
Uoo and goo- The asymptotic behavior of the solutions is described by Eq.(5.9), which 
is seen from the fact that the ratios N/Noo, Y/Yoo, and Q/Qoo approach unity, with 
Nqo, Too, Qoo dehned in (5.9). After a time reparametrization with a constant scale 
factor, the g-geometry in the far held zone is described by 

dsl = -QHR + ^ + R^dn‘^ , (7.6) 
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Figure 3. Profiles of the W2 wormhole solution obtained for the parameter value (7.5), with 
Noo, Yoo, Qoo defined by (5.9). 


where in the leading 0{B?) order the Q, N amplitudes coincide with each other and 
with the corresponding amplitude for the AdS geometry for the negative cosmological 
constant A = Kih^, 

= -^R^ + 0{R), = -^R^ + 0{R). (7.7) 

o o 

However, already in the hrst 0{R) subleading order the amplitudes Q and N are no 
longer the same and deviate from the AdS value. The f-geometry in the asymptotic 
region is expressed by Eq.(8.17) below and corresponds to a direct product x 

We did not hnd other global solutions than those of the described above types 
W1 and W2. For generic values of the parameters ki, K 2 , C 3 , C 4 , h we either do not hnd 
any solutions at all or obtain singular solutions. If the parameters are properly chosen 
and the solutions exist and extend to large values of r, then they are always found to 
be either W1 or W2. In addition, due to the symmetry (4.15), there are also solutions 
for which the two metrics are interchanged. A systematic study of the topography of 
the parameter space to identify all parameter values for which the solutions exist is a 
time consuming task that we leave for a future project. 

8 Geometry of the solutions 

Let us consider the geometry of the solutions and its global structure, in order to see 
if the wormholes are traversable or not. 

8.1 Type Wla wormholes 

Let us hrst consider solutions of the type shown in Fig.l. One can represent the 2D 
part of the g-metric as 

ds^g = —Q^dt^ + dr^ = (y—dt^ + dp^) = Q^ds‘^ . ( 8 . 1 ) 
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Here r G (—cxo, +cxo) is the proper distance, while the conformal radial coordinate 


P = 



( 8 . 2 ) 


changes within a hnite interval, p G {—poo, Poo)- The lightlike geodesics are the same 
in the dSg and ds"^ geometries. Using the Hamilton-Jacobi eqnation, the radial timelike 
geodesics followed by particles of mass p are described by 


^l] = I 

dt 


(8.3) 


where S is the particle energy. Dropping the conformal factor in (8.1) leads to the 
conformal diagram of the g-geometry in the t, p coordinates shown in Fig. 4. The central 



Figure 4. Conformal structure of the g-geometry for the Wla solution (left) and the effective 
potential in the geodesic equation (8.3) (right). 

part of the diagram, p = 0, is the position of the throat aronnd which timelike geodesics 
oscillate. As is seen in Fig.4, the effective potential has a minimnm at the throat 
position, so that the throat attracts the particles. The maximal and minimal valnes 
of the radial coordinate, ±poo, correspond to the position of the conformal timelike 
bonndary, which is at an inhnite proper distance away from the throat. Timelike 
geodesics are trapped by the conhning potential and cannot reach the bonndary, while 
nnll geodesics {p = 0) reach it in a hnite coordinate time t bnt in an inhnite affine 
time. As a resnlt, the g-geometry is geodesically complete. All this is very similar to 
the properties of the AdS geometry, apart from the fact that the bonndary consists 
now of two components, J_|_ and J_. The conclnsion is that timelike geodesics traverse 
the wormhole and oscillate aronnd the throat. 

8.2 Type Wlb wormholes 

Let ns consider the g-geometry shown in Fig. 2. The specialty now is that the amplitnde 
Q vanishes at r = r±, where r_ = —r+, so that Q is positive for r G (r_,r+) and is 
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(8.4) 


negative otherwise. The 2D part of the metric can be expressed as 


dSg = Q {—dr + dp^) = —dr^, 

where the p-coordinate is dehned by (8.2) for r G (r_,r_|_), otherwise one has 

r dr 


P = 


dr 

. Q 


and 


P = 


1-0. Q 


— if r < r_. 


This determines three coordinate regions; 

A : r G (r_,r+), 




r G (r+,oo), 
r G (-cx),r_), 


p G (—C) 0 , oo), 
P 6 (O.oo), 

P e (-00,0), 


(8,5) 


( 8 . 6 ) 


and in the last two regions p changes in the opposite directions, so that p = 0 corre¬ 
sponds to r = ±cx). In the region A one has p ± f = tan(M±) G (—cxo, cx)), so that 



Figure 5. Conformal structure of the spacetime regions defined by Eq.(8.6). 


ds^ = Q‘^{—dt^ -I- dp^) 


cos^(m+) cos^(m_) 


duj^ du- , 


(8,7) 


where u± G (—7r/2,7r/2). Therefore, the A region is conformally eqnivalent to the 
diamond in the (m+,m_) plane shown in Fig.5. The vertical symmetry axis of the 
diamond corresponds to the throat position, p = r = 0. In the B± regions one has 
either p > 0 or p < 0, hence conformal images of these regions can be obtained by 
cntting the diamond and keeping either only its right or only its left triangular part, 
as shown in Fig.5. The vertical p = 0 side of the triangles then corresponds either to 
r = CX) or to r = —oo, which is the position of the tiniehke AdS boundary. 

The null boundaries of the A, B± regions are the Killing horizons. They corre¬ 
spond to p = ±x but they can be reached by timelike geodesics in a hnite proper 
time. Specihcally, the radial timelike geodesics are described by equation (8.3) which 
can be represented in the equivalent form 


dr\^ 8"^ 

dr J 


- 1 , 


( 8 . 8 ) 
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where r is the proper time. Let us denote x = r — r+. The amplitude Q has a simple 
zero at r = r+ and close to this point one has Q = ax + 0{x^) with a constant a, in 
which case Eq.( 8 . 8 ) yields 


xdx (X dr oc (tq — r), (8.9) 

where Tq is an integration constant. This shows that starting in the A region where 
X < 0, the geodesics arrive at the boundary where x = 0 at a hnite moment of 
the proper time, r = tq. Therefore, the A region is geodesically incomplete and the 
geodesics arrive at its null boundary in a hnite proper time. The same applies to the 
regions B±. 

However, it turns out that the boundaries of the regions - the Killing horizons, 
are singular since the curvature diverges there. Introducing the orthonormal tetrad 
consisting of the vectors 

Id d Id Id 

Q dt' dr' R d’d' R sin dip ’ 

the following tetrad components of the curvature do not vanish (with ' = d/dr) 


( 8 . 10 ) 


_ Q" _ _ Q'R' 

RoiOl — -Q, -R0202 — -Ro 303 — , 


_ ^/2 

-R 2323 = -^2424 = -R 3434 = --• (8-11) 


Since Q oc x at the horizon, the components R 0202 = .R 0303 oc 1/x oc l/y'ro — r diverge, 
although this divergence is relatively mild and only leads to hnite relative deviations 
of the neighboring geodesics, and hence to hnite tidal deformations. This can be seen 
by using the equation of the geodesic deviation and integrating over the proper time 
r. However, the component i?oioi = Q”/Q diverges since Q” 7 ^ 0 when Q —)■ 0 , 
hence even the 2D metric = —Q‘^dR + is singular, even though it reduces in 
the leading order to the hat Rindler metric dSg = —a^x'^dR + dr^. Therefore, every 
radial geodesic approaching the horizon hits the curvature singularity. 

One can nevertheless try and extend the geodesics beyond the horizon in a contin¬ 
uous way, which would allow one to construct a Kruskal-type extension of the metric, 
although only within the class. Geodesics of the extended metric, after having 
crossed the horizon, enter a “T-region” where the x^ in Eq.(8.9) formally becomes 
negative, because space and time interchange their role and the metric becomes 


dSg = +Q'^{r)dt^ — dr'^ + R^{r)dVP‘. ( 8 . 12 ) 

Since r becomes the timelike coordinate, this metric describes not the static wormhole 
but rather a dynamical cosmology, a minimum of R[r) then corresponding not to the 
wormhole throat but rather to something like a cosmological bounce. A knowledge 
of such solutions would allow one to construct a maximal extension of the spacetime 
geometry in order to hnd out if the wormhole can be traversed by geodesics or not. 

However, solutions in the T-regions are presently not known. In addition, the 
maximal extension would only be continuous and not differentiable due to the singular 
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nature of the horizons (if there are solutions for which Q” vanishes at the horizon then 
their extension would be at least C^). We therefore do not pursue this line anymore 
and leave the problem of constructing a maximal extension for the Wlb wormhole 
geometry for a future project. 

8.3 Type W2 wormholes 

Let us now consider solutions of the type shown in Fig.3. The g-geometry is globally 
regular and far away from the throat is described by Eq.(5.9), so that in the leading 
order it approaches the AdS metric with the cosmological constant A = —ki6o < 0 . As 
a result, the structure of the g-geometry is essentially the same as for the type Wla 
solutions. The metric can be cast to the form (8.1), where the conformal coordinate 
p is dehned by (8.2) and changes within a hnite interval, p G (— Poo, +Poo)- This gives 



Figure 6. The structure of the g-geometry for the W2 solution shown in Fig.3 (left), and 
the corresponding effective potential in the geodesic equation (8.3) (right). 


the conformal diagram shown in Fig. 6, which is similar to that in Fig.4. The timelike 
geodesics are described by (8.3), whose potential is shown in Fig.6. The wormhole 
throat is repulsive, and in addition there is an inhnite repulsive barrier at the timelike 
boundary. As a result, particles with £ < p oscillate between the throat and the 
boundary, while those with £ > p traverse the throat and oscillate between the right 
and left boundaries, as shown in Fig.6. 

Let us now consider the f-geometry. Far away from the throat it is described by 
Eq.(5.9), so that one has asymptotically 


Y = const. xR^ + OiR), ^ ^ > 

with a constant Uoo- Let us represent the metric as 

ds'j = —q^dt^ + -I- U'^dVt^ = q^{—dt^ + dp^) + U'^dVt^ 


(8.13) 


(8.14) 


20 









with 


(8.15) 


P = 


U' 


dr. 


Jo qy 

This radial coordinate changes within a hnite range, p G {—Pooy Poo), because at large 
r one has in view of (8.13) dp oc dR/R^, therefore the integral in (8.15) converges at 
the upper limit to a hnite value poo- The radial geodesics of the f-metric obey 


dpV 2 _ T 

Tt) 


( 8 . 16 ) 


As is seen in Fig.3, the q amplitude interpolates between g(0) and goo- Therefore, 
for large enough S there are geodesics which cross the whole range of p and arrive at 
p = ±Poo in cl hnite proper time. When p —)■ ±Poo fhe geometry becomes 

ds^j = ql^{-dt^ + dp^)+ Ul^dVL^ , (8.17) 

which is completely regular. As a result, nothing prevents the f-geodesics from extend¬ 
ing beyond the values p = ±Poo- Hence, from the f-geometry viewpoint, the manifold 
corresponding to the interval p G {—poo, Poo) is geodesically incomplete, so that the f- 
geometry could be extended beyond this interval. However, as far as the g-geometry is 
concerned, the manifold is complete, because the limiting values p = ±Poo correspond 
to the AdS boundary. We therefore have a peculiar situation where the same manifold 
is complete in one geometry but is incomplete in the other. One could in principle 
try and extend the manifold by integrating the equations beyond p = ±Poo until the 
f-geometry is complete. However, the additional parts of the manifold obtained in 
this way would then be g-geodesically disconnected from the original wormhole, be¬ 
cause the latter is already g-complete. We therefore adopt the viewpoint that only the 
g-metric describes the spacetime geometry, while the f-metric should be viewed as a 
spin-2 tensor held whose geometric interpretation is possible but not necessary. 


9 Concluding remarks 

The above analysis gives strong (numerical) evidence in favor of the existence of worm- 
holes in the bigravity theory. These wormholes are very large, with the throat radius 
of the order of the inverse graviton mass, and they can be of two principal types, which 
we call W1 and W2. 

The W1 wormholes are asymptotically AdS. This feature can be understood by 
noting that the AdS space is an attractor at large r, which means the following. The 
solutions can be obtained by integrating the system of three hrst order equations (4.19) 
for N{R), Y{R), U{R). At large R the solutions approach the AdS values, so that 

Ar = A^o X (1 + z/), F = iVox(l + 0, f/ = Ai?x (1 + x), (9.1) 

where Nq = 1 — AR^/3 and the deviations ure small. Let us consider hrst 

the Wla solution shown in Fig.l. Then one has A = —0.170 and A = 0.358 (see 
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Eqs.(7.2),(7.3)). Linearizing the eqnations with respect to small z/,then gives the 
solution 

3 

v ^ ^ X ^ with s = —3, —- ± a; X i, (9.2) 

where, for the parameter values in (7.1), one finds uj = 2.068. The three different values 
of s correspond to three independent solutions, all of them approaching zero as 72 ^ cxo. 
Therefore, the stable manifold around the AdS fixed point is three-dimensional, so that 
solutions of the three-dimensional system (4.19) generically run into this hxed point, 
which is why this is an attractor. For comparison, the flat space is not an attractor 
since the stable manifold around it is only two-dimensional and the solutions miss it, 
hence they are not asymptotically flat (probably not even in exceptional cases; see 
below). 

Eq.(9.2) determines the deviation from the AdS asymptotic, 5N'^ = — Nq, 

5N'^ =+ A^/Rcos{uj\n{R) + a), (9.3) 

R 

where M, A, a are integration constants. The hrst term on the right here is the con¬ 
tribution of the massless graviton, while the second term is the effect of the scalar 
polarization of the massive graviton. The embarrassing observation is that the mas¬ 
sive contribution oscillates (this is conhrmed by the numerics) since s given by (9.2) 
has a non-vanishing imaginary part u. This indicates that the mass is imaginary. 
Indeed, the Fierz-Pauli graviton mass for fluctuations around the proportional AdS 
background is given by (in units of m) [33] 

"^FP = + y) • (9-4) 

For A = 0.358 this gives m|p = —0.37, hence the gravitons indeed behave as tachyons. 
The value of the graviton mass actually agrees with the value of uj given above, which 
can be seen by noting that the scalar graviton behaves as a scalar held. On the other 
hand, a static, spherically symmetric scalar held of mass /r on the AdS background 
decays asymptotically as 72® with 


s = 



(9,5) 


Setting here /r = 0 and choosing the minus sign yields s = —3, while setting fR = 
rripp = —0.37 yields s = —3/2 ± 2.068 x i. This reproduces precisely the values in 
(9.2). 

At the same time, one should stress that the very existence of tachyons in the 
AdS space is not necessarily a bad feature, as long as their mass squared exceeds the 
BF bound rn|p = | A (the mass for which the square root in (9.5) vanishes), in which 
case they do not produce an instability [36]. However, for the Wla solution shown in 
Fig.l one has 

3 

type Wla: mpp = —0.37 < m|p = - A = —0.12, (9.6) 
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therefore, the BF bound is violated, which implies that the solution is unstable. It 
turns out that the BF bound is violated for all Wla solutions that we could hud. 

On the other hand, for the Wlb solution shown in Fig. 2 one obtains 

3 

type Wlb: rupp = —6.01 > m^p = - A = —6.36, (9.7) 

so that the BF bound is fulhlled, therefore the tachyon instability should be absent. 
This does not immediately imply that the Wlb solutions are stable. However, since 
they do not suffer from the most dangerous instability, there is a chance that they 
could be stable, which however can only be decided after a special analysis. 

Let us hnally consider the W2 wormholes. They do not approach the proportional 
background and so it is less clear [33] how to compute the Fierz-Pauli mass. However, 
the linearization of the held equations around the asymptotic values, similar to that 
described by Eq.(9.1), gives for the deviations power law solutions with real 

powers. Therefore, there is no evidence for tachyons, so that the W2 solutions could 
perhaps be stable. It should however be again emphasized that in all cases a detailed 
stability analysis remains an open issue. 

The tachyons [30, 31] and superluminal waves [37-39] were previously detected in 
the massive gravity theory with a hxed f-metric. Their existence does not necessarily 
mean that the theory is ill-dehned but rather shows that it can have unphysical solu¬ 
tions. It seems that in the bigravity theory the situation is similar - solutions can be 
physical and unphysical [40]. The described above Wla wormholes apparently belong 
to the latter category because they show tachyons and are unstable. One should also 
say that the solutions may admit a holographic interpretation, similarly to the massive 
gravity solutions used in the holographic conductivity models [41, 42]. 

It is instructive to compare the wormholes and black holes [32]. In both cases one 
can use the Schwarzschild coordinate, = —Q‘^dt^ -|- dR^/N^ -|- R^dVL^. For black 
holes both N'^ and vanish at R = h (horizon), while for wormholes vanishes at 
R = h (throat) but does not. The bigravity black holes [32] are characterized by 

two independent values, h and a = U{h)/h, and they can be obtained by integrating 
Eqs.(4.19) for N{R), Y{R), U{R) with the boundary condition N{h) = Y{h) = 0. The 
equation Q' = FQ (Eq.(4.20)) then insures that Q{h) = 0, since one generically has 
at the horizon 2F = 1/{R — h) + 0{1). Now, the wormholes are actually the same 
solutions but obtained for special values of a (given by Eq.(6.7)) for which the pole of F 
is canceled, and so the equation Q' = FQ ensures that Q is hnite as R = h. From this 
viewpoint, wormholes can be viewed as the special case of black holes corresponding 
to the hne-tuned a. 

The bigravity black holes generically approach the AdS space [32], but in ex¬ 
ceptional cases, for specially adjusted values of a (and for h > 0.86), they can be 
asymptotically flat [43]. For the wormholes the value of a is already hxed by the con¬ 
dition of having a regular throat, so that one cannot further adjust it to fulhll the 
asymptotic hatness condition as well. Therefore, asymptotically hat wormholes are 
unlikely to exist. 

The symmetric wormholes exist only in the bigravity and not in the massive 
gravity theory with a hat f-metric. Indeed, the hat f-metric requires that E = 1, 
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which is not compatible with the boundary condition expressed by (6.2). However, we 
have checked that in the massive gravity limit there are non-symmetric under r —)■ —r 
wormhole-type solutions for which R develops a minimum, and even inhnitely many 
minima. 
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